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PREFACE 


The  purpose  of  this  Memorandum  is  to  present  a  comprehensive  and 
accurate  table  of  values  for  the  similar  solutions  of  the  laminar 
boundary-layer  equations  for  a  homogeneous  gas.  These  values  may  be 
used  to  determine  skin  friction,  heat  transfer,  and  appropriate  bound¬ 
ary-layer  thickness  parameters  for  a  wide  'ariety  of  physical  situa¬ 
tions.  They  are  applicable  to  such  problems  as  stagnation-point  flo..s, 
flow  over  a  wedge  or  a  flat  plate,  and  hypersonic  viscous  interaction, 
and  may  also  be  used  as  good  approximation  for  the  determination  of 
nonsimilar  flows.  Well  over  1800  solutions  are  tabulated,  of  which 
about  five -sixths  have  not  been  published  before.  The  study  is  part 
of  RAND's  work  for  ARPA  on  reentry  aerodynamics.  C.  Forbes  Dewey,  Jr. 
is  an  Assistant  Professor  in  the  Department  of  Aerospace  Engineering 
Sciences  at  the  University  of  Colorado  and  a  consultant  to  The  RAND 


Co rp oration. 
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SUMMARY 


The  purpose  of  this  Memorandum  Is  to  present  a  comprehensive  and 
accurate  table  of  values  for  the  similar  solutions  of  the  laminar 
boundary- layer  equations  for  a  homogeneous  gas.  These  values  may  be 
used  to  determine  skin-friction,  heat-transfer,  and  appropriate  bound¬ 
ary-layer-thickness  parameters  for  a  wide  variety  of  physical  situa¬ 
tions.  They  are  applicable  to  such  problems  as  stagnation-point  flows, 
flow  over  a  wedge  or  a  flat  plate,  and  hypersonic  viscous  interaction, 
and  they  may  also  be  used  as  good  approximations  for  the  determination 
of  nonsimilar  flows. 

The  general  boundary- layer  equations  for  a  thermally  and  calori- 
cally  perfect  binary  gas  mixture  are  derived  together  with  a  discus¬ 
sion  of  the  numerical  integration  procedure  used  to  obtain  solutions 
to  the  equations,  the  concept  of  local  similarity,  and  solutions  for 
large  values  of  the  pressure-gradient  parameter. 

The  tabulated  solutions  include  effects  of  pressure  gradient 

2 

(0  <■  3  ^  5  and  P  **  ®) ,  Mach  number  (0  5  U  /2H  S  1)  temperature- 

00  0 

viscosity  law  (0.5  <  *  <  1  or  0  <  s  <  0,2),  leading-edge  sweep 

,0. 1  ^  t  ^  1)  wall  temperature  (0  ^  t  ^2),  Prandtl  number 
s  w 

2  2 

(0.5  <  Pr  ^  1) ,  local  streamwise  velocity  (0  ^  u  /u  -  1),  and 

e  °° 

mas;-*  transfer  at  the  surface  (-1.6  ^  f  ^  0). 

w 
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SYMBOLS 

A.  ...  A,  *  constants  defined  in  Eqs.  (105),  (115),  and  (120) 

1  D 

a  =  function  defined  in  Eq.  (109);  also  function  defined  in 
Eq.  (154) 

B  =  function  defined  in  Eq.  (86) 

C  ■  Chapman-Rubesin  constant 
C(a)  ■  limit  of  g(a,t)  as  t  -•  * 

o 

C-  *  local  skin-friction  coefficient;  C,  *»  2t/p  U 

f  f  CO  OD 

C.  *  local  heat  transfer  coefficient;  q  /p  U  (H  -  H  ) 
h  n*  ®  ®  aw  w 

th 

c^  =  mass  concentration  of  i  component 
0^  =  specific  heat  of  the  iC^  component 

=  binary  diffusion  coefficient 

E  =»  function  as  defined  in  Eq.  (102) 

£  =  numerical  difference  function  defined  in  Eq.  (82) 

«  =  function  a3  defined  in  Eq.  (86) 

=  similarity  function  defined  in  Eq.  (23) 

f  ‘  =  transformed  fluid  velocity  as  defined  in  Eq.  (12) 

fw  ■»  wall  constant  defined  by  Eq.  (22) 

3  =  function  defined  in  Eq.  (83) 

G  =*  dimensionless  enthalpy  function,  (H/He) 

g  ■  transverse-velocity  function,  w/w^ 

S  =  numerical  difference  function  defined  in  Eq.  (84) 

H  =»  total  enthalpy  of  mixture  ■  h  +  (l/2)(u  +  w  ) 
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1^3) 

Jn(t) 


j 

k 

Le 

M 

in 

Pr 

P 

4 
R 

r 

5 
Sc 

s 


s 

o 


numerical  difference  function  defined  in  Eq.  (85) 
chemical  enthalpy  of  i1"^  component 


integrals  defined  by  Eqs.  (74)  to  (78) 


th 

function  defined  as  the  n  11  integral  of  the  error  func¬ 
tion 

geometrical  index  in  boundary  layer  equations ;  also 

Reynolds  analogy  ratio  2C  /C- 

n  f 

thermal  conductivity  of  the  mixture 

Lewis  number  p  D,_  c  /k 
12  p 

Mach  number 

constant  exponent  defined  in  Eq.  (34) 

Prandtl  number,  c  n/k 
P 

fluid  pressure 
heat  flow 
gas  constant 

radial  coordinate  define  in  Fig.  1 
Sutherland  constant 
Schmidt  number,  p./p 

S/T 

o 

transformed  velocity  function  defined  in  Eq.  (150) 
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T  ■  temperature  of  the  fluid 

m  free-stream  stagnation  temperature 

t  ■  transformed  similarity  variable  defined  in  Eq.  (150);  trans¬ 
formed  variable  defined  by  Eq.  (108) 

t  -  dimensionless  adiabatic  wall  temperature,  t  ■  t  for  8 '(0)  ■  0 
aw  fw 

tr  ■  Eckert  reference  temperature  defined  in  Eq.  (53c) 

tg  «  dimensionless  sweep  parameter  defined  by  bq,  (32) 

-  1  -  (U2/2H  )  sin2  A 

oo  e 

t  •  dimensionless  wall  enthalpy  ratio  defined  by  Eq.  (31)  ■  T^/T q 

U  ■  free-stream  velocity 
00 

u  *  flow  velocity  in  the  x-direction 
v  ■  flow  velocity  in  the  y-direction 
W  ■  constant  defined  in  Eq.  (141) 

w  «  flow  velocity  in  the  z-directien  (transverse  velocity) 
x  -  coordinate  in  direction  of  flow 
y  ■  coordinate  normal  to  the  surface 
Z^  ■  dimensionless  concentration  function  defined  in  Eq.  (15) 

z  **  dimensionless  function,  (1  -  Z^)  ;  also,  coordinate  transverse 
to  the  flow  direction 

8  ■  pressure  gradient  parameter  defined  in  Lq.  (34b) 

Y  ■  adiabatic  constant  ■  c  /c 

P  v 

* 

6  -  boundary  layer  displacement  thickness  defined  by  Eq.  (80) 

e  *  asymptotic  parameter  defined  in  Eq.  (98) 

7)  -  similarity  variable  defined  in  Eq.  (11) 

®  ■  boundary- layer  momentum  thickness  defined  in  Eq.  (81) 


xiv 


8  «  dimensionless  enthalpy  function,  (H  -  H  )/(H  -  H  ) 

w  e  w 

h  -  trial  function  for  solution  of  first-order  equation  de¬ 
fined  by  Eq.  (122) 

A  -  sweep  angle 

X  ■  dimensionless  density-viscosity  product,  (l/C)(pM>/p  |A  ) 

6  C 

pi  ■  viscosity 

§  »  transformed  x-coordinate  defined  in  Eq*  (10) 

P  "  fluid  density 

2 

c  ■  hypersonic  parameter  (U  /2H  ) 

GO  0 

2  2 

a.  ■  modified  hypersonic  -*aran:ater,  (U  /2H  )  •  (u  /u  ) 

1  GO  0  0  GO 

a„,  x  modified  hypersonic  parameter, 

-  (V^/2H  )[{u  /u  )2  cos2  A  +  sin2A] 

-  e  e  * 

t  -  total  skin  friction 

x  skin  friction  in  spanwise  direction 

t  ■  skin  friction  in  x-direction 
x 

$  x  function  defined  in  Eq.  (105) 

<  ■  variable  defined  by  Eq.  (141) 

X(a)  x  function  defined  in  Eq.  (112) 

XQ  "  constant  defined  in  Eq.  (141) 

K*.y)  ■  stream  function  defined  by  Eq.  (9) 

♦l  ...  ■  functions  defined  in  Eqs.  (117),  (118),  and  (121) 

ts  x  exponent  in  the  viscosity-temperature  law,  pi  ~  Tu' 
a>r  x  constant  defined  <n  Eq.  (53b) 


<ffwnrn*wfw*ityff>wni  smwmngwtrym*r» 


-XV- 


Subscripts 

(  )  *  function  at  the  edge  of  the  boundary  layer 

(  )  ■  function  at  the  wall 

w 

(  )  =  function  evaluated  in  the  free  stream 
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1.  INTRODUCTION 


Daring  the  last  decade,  many  similar  solutions  to  the  laminar 
boundary- layer  equations  have  been  obtained.  Most  of  these  results 
have  been  for  specific  exact  conditions,  such  as  stagnation-point 
flow  or  flat-plate  flow,  and  are  based  on  the  assumptions  of  a  linear 
temperature-viscosity  law  and  a  Prandtl  number,  Pr,  of  unity.  It  is 
well  recognized,  however,  (cf.  the  citations  of  Refs.  1  to  5)  that  an 
accurate  estimate  of  skin  friction,  heat  transfer,  and  boundary- layer 
thickness  in  compressible  flow  usually  requires  the  use  of  a  realis¬ 
tic  temperature-viscosity  law  and  retention  of  the  dissipation  term 
which  appears  in  the  energy  equation. 

The  purpose  of  this  Memorandum  is  to  present  comprehensive,  sys¬ 
tematic,  and  accurate  tables  of  solutions  to  the  laminar  boundary- 
layer  similarity  equations  for  a  perfect  homogeneous  gas.  These  solu¬ 
tions  are  applicable  to  the  classical  cases  wherein  the  requirements 
for  similarity  are  satisfied  exactly.  The  solutions  may  be  used  also 
in  applying  local  similarity  methods  to  cases  that  do  not  meet  the  ex¬ 
act  requirements  for  similar  flows.  In  particular,  the  sensitivity  of 
the  numerical  values  of  the  wall  derivatives  (which  govern  heat  trans¬ 
fer  and  skin  friction)  and  the  values  of  the  integral  thickness  parame¬ 
ters  to  changes  in  the  similarity  variables  may  be  accurately  determined. 
For  many  physical  situations,  relaxing  one  or  more  of  the  conditions 
required  for  exact  similarity  will  not  lead  to  large  errors  when  the 
local  similar  solutions  are  used  to  estimate  boundary- layer  properties. 

Solutions  are  included  to  illustrate  the  effects  of  ieading-edge 
sweep,  mass  transfer  at  the  surface,  pressure  gradient,  wall  temperature. 
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free-stream  Mach  number,  local  external  flow  velocity,  the  Prandtl  num¬ 
ber,  and  the  viscosity-temperature  law  of  the  fluid.  The  tables  in¬ 
clude  most  of  the  numerical  solutions  obtained  by  previous  authors 
(recomputed  to  the  accuracy  of  the  present  program)  and  also  new  so¬ 
lutions.  r"hese  tables  do  not  include  solutions  with  3  <  0  (decelerat¬ 
ing  flows)  or  f  >  0  (suction). 

Section  11  presents  the  similarity  equations  in  their  general 
form,  including  the  equation  for  species  concentration.  Transforma¬ 
tions  applicable  to  compressible  flows  are  given,  as  are  equations 
for  the  computation  of  skin  friction,  heat  transfer,  and  integral 
boundary- layer  thicknesses.  Classes  of  reduced  equations  are  also 
discussed;  these  classes  represent  cases  where  one  or  more  of  the  pa¬ 
rameters  Pr,  is,  (3,  t  ,  t  ,  IT/2K  and  u^/u^  are  either  zero  or  unity. 

Section  III  gives  a  concise  guide  to  the  solutions  of  this  Memo¬ 
randum.  Section  IV  presents  the  numerical  method  of  computation  used 
in  this  program.  The  concept  of  local  similarity  is  discussed  in  Sec¬ 
tion  V,  where  criteria  are  presented  for  judging  the  applicability  of 
similar  solutions  in  situations  where  the  similarity  parameters  vary. 

The  analysis  is  extended  to  large  positive  values  of  the  pressure-gra¬ 
dient  parameter,  3,  in  Section  VI. 

Our  experience  indicates  that  values  of  the  skin-friction  and  heat- 
transfer  derivatives  and  the  boundary- layer  thickness  integrals  may  be 
estimated  to  high  accuracy  for  intermediate  values  of  the  similarity 
parameters  by  careful  cross-plotting  of  the  present  exact  solutions. 

The  accuracy  varies  witn  the  particular  quantity  and  the  specific  val¬ 
ues  of  the  similarity  parameters,  but  it  is  generally  between  0.5  per¬ 
cent  and  2.0  percent.  In  particular,  use  of  the  solutions  for  f3  -*  ® 
allows  estimation  of  all  quantities  for  all  positive  3. 
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II.  SIMILARITY  EQUATIONS  FttR  THE  LAMINAR  BOUNDARY  LAYER 

The  purposes  of  this  section  are  three.  First,  we  shall  display 
the  laminar- boundary- layer  equations  for  a  binary  mixture  of  perfect 
gases  in  their  most  general  form,  lis“fng  the  general  conditions  that 
lead  to  similar  solutions.  Second,  the  specific  equations  treated  here 
are  deduced  from  the  general  equations,  and  relations  are  given  for  com¬ 
puting  heat  transfer,  skin  friction,  and  integral-thickness  parameters. 
Finally,  special  reduced  forms  of  the  general  equations  (Pr  =  1,  etc.) 
are  discussed. 


GENERAL  EQUATIONS 

The  equations  that  describe  the  conservation  of  mass,  momentum, 
and  energy  in  a  lamin'. r  boundary  layer  are  well  known.  They  can  be  ex¬ 
pressed  as  follows: 


—  (purJ)  +  —  (pvrJ)  =  0 

dx  dy 


(1) 


Du  du 

pu  —  +  pv  — 
dx  dy 


dw  dw 

pu  —  +  pv  —  “ 

dx  dy 


(2) 


(3) 


dp 

—  =  0 
dy 


(4) 
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9c. 
i  , 

pu  - +  pv 

9x 

9c.  Id/.  9c.  \ 

-1  -  —  —  (t'p  »12  — ) 

9y  rJ  Sy  \  dy  / 

(6) 

j  =  o 

two-dimensional  flow 

(7a) 

j  -  1 

axi symmetric  flow 

(7b) 

The  coordinates  x,  y,  and  r  are  defined  in  Fig.  1.  Other  symbols  are 
defined  as  follows: 

*  mass  fraction  of  i*"^  constituent 

Dj^  *  binary-diffusion  coefficient 

H  =  total  enthalpy  of  mixture  =  h  +  (l/2)(u^  +  w^) 

til 

tu  =  chemical  enthalpy  of  i  constituent 

k  =  thermal  conductivity  of  fluid 
Le  =  Lewis  number,  P  cp/k 

Pr  =  Prandtl  number,  c  p./k 

P 

p  «=  local  fluid  pressure 
Sc  3  Schmidt  number,  \i/p 

u,v  =  flow  velocities  in  the  x,  y  directions,  respectively 
w  =  flow  velocity  in  direction  normal  to  x,  y  plane 
p,  =  viscosity  coefficient  of  fluid 
p  =  local  fluid  density 
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These  equations  neglect  thermal  diffusion  and  diffusion-thermo  effects. 

If  the  fluid  is  a  mixture  of  thermally  perfect  gases,  Eq.  (4)  in¬ 
tegrates  to 


p(x,  y)  -  Pe(x)  =  PeRTe 

(8) 

A  stream  function  Y(x,y)  is  introduced  that  satisfies 

equation: 

the  continuity 

dY/dx  -  -  pvr^ 

(9a) 

SY  /  by  *  +  pur^ 

(9b) 

and  a  transformation  for  the  independent  variables  is  introduced: 


I 


r 


Cp  h*  <*’ 
e  e  2  k 


r.2> 


dx 


(10) 


T] 


(ID 


where  subscript  e  denotes  any  property  at  the  edge  of  the  boundary  lay¬ 
er  and  r^  denotes  a  characteristic  radius.  The  quantity  C,  by  defi¬ 
nition,  is  a  combination  of  physical  properties  that  are  functions  of 
x  only. 

o 

In  reducing  the  partial  differential  equations  in  1}  and  |  to  or¬ 
dinary  differential  equations  in  the  single  variable  T|,  all  derivatives 
with  respect  to  §  within  the  boundary  layer  are  neglected.  The  fluid 
velocity,  the  enthalpy,  and  the  concentration  are  then  defined  in  terms 
of  the  transformed  similarity  variable  T|: 
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-  df  (Tjj  / dT|  =  f  ;(T|) 

(12) 

w/w 

e 

-  gOl) 

(13) 

H/H 

e 

-  G  (Tj) 

(14) 

c.  /c. 
i  iw 

-  Zt(Tl) 

(15) 

(Tl*a  dependent  variable  f  appearing  in  Eq.  (12)  is  simply  equal  to 
/31b  times  the  stream  function  i.)  Then  the  boundary- layer  equations, 
(2),  (3),  (5),  and  (6)  for  a  two-component  boundary  layer  may  be 
written  as  a  set  of  coupled  ordinary  nonlinear  differential  equations. 
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where 

A  ™  leading-edge  sweep  angle 

X  «  (l/C)(p$i/p  n  ) 
e  e 

u^  -  free-stream  flow  velocity  in  the  x-direction 
The  above  equations  are  applicable  to  the  two-component  boundary 
layer  of  an  axi symmetric  body  or  a  swept  two-dimensional  surface.  The 
boundary  conditions  for  these  equations  are  given  at  the  surface  of  the 
body  by  the  no-slip  condition  and  the  requirement  of  similarity  for 
the  wail  concentration  and  enthalpy: 


f  '(0) 

=  0 

(20a) 

8(0) 

a  0 

(20b) 

G(0) 

-  G 

w 

(20c) 

Z^O) 

-  1 

(20d) 

At  the  outer  edge  of  the  boundary  layer,  the  velocity  and  enthalpy  must 
match  the  values  in  the  free  stream.  The  concentration  of  any  material 
injected  at  the  surface  is  assumed  to  be  zero  in  the  external  flow. 


f'W 

-  1 

(21a) 

g(°>) 

=  1 

(21b) 

G(°0 

=  i 

(21c) 

Zx(®) 

=  0 

(2  Id) 
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The  final  boundary  addition  is  obtained  by  considering  the  2-parcicles 
at  the  wall.  This  yields  the  well  known  Ecker t-Schneider  condition 
represented  by  the  equation: 


f  -  f(0) 

w 


Z'(0) 


(22a) 


From  the  continuity  equation,  the  injection  parameter  is  related  to 
the  similarity  coordinates  by  the  equation 


f 

w 


f 


w 


(pv)  fit 
w 


Cp  p.  u 
e  e  e 


all  0 


du 


dx 


w 


21 

CP  i  u  r.  J 
e  e  e  k 


-1/2 


>  (22b) 


P  *  0 


CONDI  HONS  FOR  SIMILARITY 

In  order  that  Eqs.  (16)  through  (19)  be  similar,  severe  restric¬ 
tions  must  be  placed  on  some  of  the  terms.  These  are: 


2§  du  .  „  P  >. 

-7f(£  -f)  '  Fi™ 

u  ds  P 


(23) 


XOl) 


l2/rk  =  F2(T1) 


(24) 

(25) 


either  Pr  =  1 


.2,  2 


(26) 


or  u  /u  =  const,  and  Pr  =  ?r(l) 

e  CD  x 


1 
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Le  =  Le("n) 


(27) 


2  2 

G  ,  c,  ,  and  r  /r,  are  constants 
w  lw  w  k 


(28) 


f  =  f(0)  =  const, 

w 


(29) 


Equations  (23)  and  (25)  represent  the  most  severe  restrictions  on  the 
system  because  the  external  velocity  distribution  and  the  body  surface 
must  satisfy  the  requirements  shown.  It  may  be  verified  that  Eq.  (23) 
has  th';  following  form: 


Fj(Tl) 


21  du  T  (  ,  1  r 

- -■  ~  t  \f'2 - f(  1  -  t 

u  d|  T  6  l  t  L  w 


)6  -  (1  -  t8)g"  +  t 


(30) 


where 

9  *  (H  -  H  }/(H  -  H  ) 

w  e  w 

Tq  =  free-streara  stagnation  '  erature 

and 

t  =  H  /H  -  T  /T  (31) 

w  w  e  w  o 

(adiabatic  and  calorically  perfect  i.  wiscid  flow)  and 


,  .  Y  -  1  -.2  2  . 

1  +  — r—  M  cos  A 
2  00 


1  + 


Y 


-in2 


(32) 


It  should  be  noted  that  either  1-0  and  0^t  ^  1,  or  j  =  1  and 
t  =1.  If  t  and  t  are  considered  to  be  free  parameters,  Eq.  (30) 

S  S  V 


requires  that 
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2?  du  T 

- £-a  t 

u  d|  T  8 
e  e 


♦1(T1) 


(33) 


If 


we  assume  Che  relationship  between  u 


e 


and  %  to  be  given  *-  r 


u 

e 


/t"7t  ira 

e  o 


(34a) 


Chen 


2§  du  /T  \ 

-rf\f  k  -  2"  *  6 

u  d§  '  e ' 
e 

where  P  is  the  edified  Falkner-Skan  parameter. 

Using  Eq.  (34)  to  define  the  parameter  (3,  the  equations  assume 
the  following  similarity  form: 
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OS) 


2 

In  Eq.  (37.',  <7  =  (U^/2He)  and  13  given  the  designation  of  hypersonic 
parameter.  Two  modified  hypersonic  parameters  can  be  defined,  as 
follows 


al 


(39a) 


and 


a 


2 


cos  A  +  sin  A 


. 


(39b) 


The  parameter  includes  the  shock-angle  effect  on  the  energy  equation 
it  is  a  times  the  ratio  of  the  velocity  at  the  edge  of  the  boundary 
layer  to  the  free-stream  velocity.  The  parameter  a ^  is  a  generaliza¬ 
tion  of  including  the  effect  of  sweep  angle.  When  a  =  0,  the  ratio 
of  kinetic  energy  in  the  flow  to  the  stagnation  enthalpy  of  the  free 
stream  is  zero;  as  a  result  the  dissipation  terra  in  the  energy  equation 
is  neglected.  As  a  increases,  the  dissipation  term  plays  an  ever-in¬ 
creasing  role  in  the  energy  equation. 

EQUATIONS  FOR  SOLUTION 

Subject  to  simplification  regarding  the  flows  to  be  considered, 
Eqs.  (35)  to  (38)  are  those  equations  for  which  solutions  are  presented 
in  this  Chapter.  These  simplifications  are: 
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1.  The  Chapman-Rubesir.  constant  is  defined  as 


C  -  p  P  /p  P 
w  w  e  e 


so  that  the  quantities  £  and  X  become 


rr.^  p  u  dx 
k  w  w  e 


X  *  pp/p  p 
w  w 


2.  Curvature  effects  are  neglected: 


(r/rk)‘  =  1 


3.  The  product  (Le  -  1)  •  (hj  -  h^)  is  assumed  to  be  zero,  and 
thermal  diffusion  is  neglected.  A  heuristic  justification  of  this 
simplification  for  equilibrium  air  is  >iven  by  Beckwith  and  Cohen, 
Appendix  A  of  Ref.  6. 

The  Lewis  number  has  been  taken  to  be  unity;  i.e. ,  Pr  *  Sc.  As 
a  final  step,  the  concentration  variable  Z^  is  transformed  to  the  de¬ 
pendent  variable 


Z  =*  1  -  Z. 


c  -  c 
lw  1 


Substitution  of  9  for  G  and  z  for  Z^  yields  the  equations  that 
were  solved  exactly  to  prepare  the  tables  of  solutions: 


(Xf")  '  +  ff" 


(1  -  tw)0  -  (1  -  ts)g‘+  twJ>  (40) 
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(XgV  '  fg'  -  0  (41) 


This  is  a  9t^1-order  nonlinear  set  of  ordinary  differential  equations. 

Nine  boundary  conditions  are  required,  five  at  the  wall  and  four  at 
the  edge  of  the  boundary  layer.  The  boundary  conditions  to  be  satis¬ 
fied  are: 

At  the  wall  when  "H  =  0: 

£*(0)  a  g(0)  «■  6(0)  ■  z(0)  =  0  (44) 

and  the  Eckert-Schneider  condition  which  relates  the  conservation  of 
2~particles  at  the  wall: 

£ (0)  - - — -  z  '(0)  (45) 

S'd  -  clv) 


In  obtaining  Eq.  (45)  from  Eq.  (22a),  we  have  made  use  of  the  fact 

that  (X)  *  1.  At  the  outer  edge  of  the  boundary  layer,  the  condi- 

w 

tions  are: 


f'{\)  =  g(y  =  9(y  =  z(7le)  =  1  (46) 

The  value  of  T|  is  chosen  to  be  large  enough  to  insure  that  the  bound¬ 
ary  conditions  given  by  dq.  (46)  are  satisfied  asymptotically  to  a  high 
degree  of  accuracy  (see  Section  IV). 
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The  similarity  conditions  are  satisfied  by: 


3  =  constant 

(47) 

X  =  pp/p  P  “  X(T1) 

w  w 

(48) 

either  Pr  =  1  | 

(49) 

or  (u  /u  )2 
e  00 

=  constant  and  Pr  =  Pr  (TJ)  j 

Sc  = 

3c(T|)  [  =  Pr(T])  for  Le  =  1] 

(50) 

t  and  c.  are  constants 

W  -W 

(51) 

f(0)  =  constant  =  f 

w 

(52) 

The  viscosity-temperature  relationship  may  be  characterized  by 
the  power-law  expression:  p  =  AT^.  This  yields  X  ~  T111  A  value 
(u  ■  0. 7  corresponds  to  conventional  wind-tunnel  conditions,  ‘.’hile 
a)  b  0.5  represents  conditions  encountered  in  hypersonic  flight. 

In  terms  of  the  similarity  variables,  X  may  be  written  in  the 

form 


X 


—  {(1  -  tw)6  -  (1  -  tg)g2  -  (cos2  A)f /2  +  tj 


If  a)  ■  1,  then  X  =  1  and  pp  =  constant.  In  the  literature,  this  lat¬ 
ter  assumption  has  often  been  made  in  order  to  simplify  the  boundary- 
layer  equations. 

(4) 

The  authors  have  demonstrated  in  a  previous  paper  that  bound¬ 
ary-layer  characteristics  calculated  by  means  of  a  Sutherland  viscosity 
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law  can  be  approximated  almost  exactly  using  the  power- law  relation 
p  ~  T  ,  provided  the  empirical  exponent  <jl'  is  suitably  chosen.  Suther¬ 
land's  law  may  be  written 


P 


o 


3/2 


s  +  1 


s  +  t 


(53a) 


where  t  =  t/T  and  s  =  S/T  .  The  Sutherland  constant,  S,  is  a  charac- 
oo 

teristic  temperature  for  the  gas,  and  p^  is  the  viscosity  evaluated  at 

the  stagnation  temperature  Tq.  The  empirical  equations  for  calculating 

a’  are 
r 


0) 

r 


3 


£n 


t  +  s 
w 

t  +  s 
r 


2  in  t 

r 


(53b) 


t  -  0.5  (t  +  t  )  +  0.22  (1  -  t  ) 
r  e  w  w 


(53c) 


The  quantity  t  is  sometimes  referred  to  as  the  Eckert  reference  tem¬ 
perature. 

Solution  of  Eqs.  (40)  to  (43)  subject  to  the  bour. ’ary  conditions 
of  Eqs.  (44)  to  (46)  requires  the  specification  of  eight  independent  pa¬ 
rameters.  These  are 

f  f(0)  =  blowing  parameter 


Pr 

t 

s 

t 

w 


Prandtl  number,  taken  to  be  constant  *  Sc 
sweep  parameter 

normalized  wall  temperature 

local  strecawise  velocity  ratio 
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3  =  pressure-gradient  parameter 

2 

a  =  U  /2H  =  hypersonic  parameter 

CD  6 

it  or  s  =  temperature-viscosity  law 

All  eight  must  be  independent  of  the  streamwise  coordinate  §  for  simi- 

2 

larity  to  hold  exactly.  If  Pr  =  1,  the  parameters  o  and  (u^/uj  are 
not  required. 

It  should  be  emphasized  that  thermal  diffusion  has  been  neglected 
and  the  product  (Le  -  l)(h^  -  h^)  has  been  assumed  zero.  In  this  ap¬ 
proximation,  solution  of  the  three  equations  for  f,  0,  and  g  are  in¬ 
dependent  of  the  diffusion  equation  and,  consequently,  are  independent 
of  the  value  of  the  Schmidt  number. 

SPECIAL  CLASSES  OF  REDUCED  EQUATIONS 

Two-Dimensional  Flow 
2i  2\ 

If  r  J/r k  =  1,  the  ordinary  two-dimensional  boundary- layer  equa¬ 
tions  result: 
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(Z  f 

X  ~M  +  fz'  -  0  (57) 

Equations  (54)  to  (57)  are  valid  for  axisymmetric  flow  if  A  is  set 
equal  to  zero  and  j  is  unity. 


A  -  0 

When  the  free-stream-f low  direction  coincides  with  the  x-direction, 

then  A  ■  0  and  t  -  1.  For  this  case,  the  boundary- layer  equations  re- 
s 

duce  to: 


(Xf') #  +  ff" 

*  e[f'2 

-  (1  -  t  )0  +  t  ] 
w  wJ 

(58) 

1  e'\' 

( /i  i 

^  2aX  /J  \  Y 

1  -  +  f9'  - 

i  (~  ■  1 

) -  (-§  Jf'  f  '  J 

\  Pr/ 

j  \Pr 

'  a  -  v  1 

1  )  (*1  -  h2>Zl'l 

H  LSc  \Le  /  \ 


(59) 


(x^  ( +  £z‘ 


(60) 


No  Mass  Transfer 

In  the  special  case  where  no  mass  transfer  takes  place  in  the 
boundary  layer,  Z^  =  0  and  c^  =  0,  and  the  boundary  conditions  (in- 

•ff 

eluding  the  Eckert-Schneider  condition'  ure  ignored: 

*Although  the  solution  for  the  concentration  ratio  is  Z^  *  0, 
the  solution  of  Eq.  (43)  for  the  normalized  concentration  z  has  a  non¬ 
vanishing  solution  and  a  nonzero  gradient  z'  at  the  wall. 
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(Xf")  '  +  ff" 


<  eV 

x  —  +  fe' 


Bf£'2  -  (1  -  tu)0  -  tu] 


2a  X 


(1  -  t  ) 
w 


f"f ' 


(61) 


(62) 


Equations  (61)  and  (62)  are  the  ordinary  two-dimensional  Prandtl  bound¬ 
ary-layer  equation. 


U  ~  T 

If  the  viscosity  is  directly  proportional  to  the  temperature, 
>  =  1,  and  Eqs.  (61)  and  (62)  reduce  to: 


f  +  ff' 


= 


a  - 


t  )e 

W 


-«} 


(63) 


(64) 


Pr  -  1 

Equations  (63)  and  (64)  reduce  to  a  form  in  which  the  dissipation 
term  of  the  energy  equation  is  zero: 


f*  +  ff"  =  0  [f  '2  -  (1  -  tw)6  -  tw  j  (65) 

0"  +  fe'  -  0  (66) 


3-0 

Finally,  if  the  pressure-gradient  parameter  is  zero,  the  momentum 
equation  is  uncoupled  from  the  energy  equation;  Eq.  (65)  becomes 


20- 


f"'  +  ££"  =  0  (67) 

and  by  inspection,  the  solution  of  the  energy  equatio-  is 

6  =  £'  (68) 

Computation  of  Boundary-Layer  Properties 

Heat-transfer  rates  and  skin  friction  are  related  by  the  similar¬ 
ity  transformations  to  the  derivatives  6/(0),  ^'(0),  and  g/(0).  The 
heat  transfer  from  the  stream  to  the  wall  is  given  by 

^  =  +kv(9I/3y)v  (69a) 


After  the  proper  substitutions  have  been  '*ede,  the  expression  for  the 
heat  transfer  in  similarity  coordinates  is 


k  p  u  H  (1  -  t  )r.J 
w  w  e  e  w  k  n 

%  m  *  - E~7!r - 6  (0) 

p» 


(69b) 


and  the  local  heat  transfer  coefficient  is 


C.  -  q  /[p  U  (H  -  H  )] 

h  vtf  oo  oo  aw  w 


(70) 


The  skin  friction  for  the  x-direction  can  be  described  in  similarity 
coordinates  as  follows: 


p 

u  i\0) 


w  m 


The  spanwise  coefficient  is 


r.^  p  u  w 
k  wee 

8  #(0) 


(71) 


(72) 
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The  total  component  of  skin  friction  is  the  vector  sum  of  the  two  com¬ 
ponents  t  and  Tz : 


rp  p  u 
k  w  e 


w 


lifwY 


U2  + 

e 


[g'«>>r 


1/2 


(73) 


and  the  local  skin-friction  coefficient  is  defined  by 


C 


f 


2i/p  UA 


(74) 


Several  integral  relations  have  also  been  tabulated.  They  are 


11  -  (1  -  t8)Cl1(l)  -  Ii(2)1  -  (1  -  tw)I1(3)  +IX(2)  (75) 

00 

12  -  J  f'(l  -  f#)  dT)  (76) 


where 


I^D 

Ix(2) 

Ix(3) 


f  (1  -  g2)  dT] 


o 


o 


(1  “  f'2)  dT] 


(1  -  0)  dT) 


(77) 


(78) 


(79) 


These  integrals  can  be  related  to  the  boundary- layer-thickness  parame¬ 
ters  expressed  in  similarity  coordinates. 

Displacement  thickness: 
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r00  /  pu  \  ^2?  /T  V 

/Te\ 

1  1 1  -  j  dy  =  .  1  j 

‘°  ’  W  Wk'V 

''-(th 

L  O 

Mementum  thickness : 


(80) 
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P  u  r.^ 
e  e  k 


(81) 


RANGE  OF  SOLUTIONS  AND  PARAMETERS 


The  parameters  for  the  systems  of  equations  under  consideration 
here  are 

pressure-gradient  parameter  0  5  p  £  5  and  8  *  <*> 

Prandtl  number  Pr  =  0.5,  0.7,  1.0 

Schmidt  number  Sc  =  Pr  (i.e. ,  Le  =■  1) 

temperature-viscosity- law  parameter  id  =  0.5,  0.7,  1.0  or 

.01  5  s  5  .3 

sweep-angle  parameter  0  £  t  5  ? 

s 


mass-transfer  parameter 
wall  temperature 
shock-angle  parameter 
hypei conic  parameter 


-0.6  5  f  ^0 
w 

0  5  t  5  1.2 
w 

(u  /u  )2  -  0,  0.5,  1.0 
0  00 

a  ■  IT/2H  r.  0,  0.5,  1.0 

oo  e 


The  quantities  tabulated  are 

f^  =  normalized  chordwise  velocity  gradient  at  wall,  [f7/(0)"J 

0'  =  normalized  total  enthrlpy  gradient  at  wall,  [07(O)] 

w 

g7  »  normalized  transverse  velocity  gradient  at  wall,  Tg,(0)j 
w  ° 
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plus  the  integrals  1^  1 2%  1^1),  1^2),  1^3)  as  defined  by  Eqs.  (75) 
to  (79). 
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III ,  NUMERICAL  RESULTS 


This  section  provides  a  guide  to  the  present  solutions  of  the 
laminar  boundary- layer  equations.  The  available  solutions  are  shown 
schematically  in  Keys  1  to  6.  Each  key  shows  all  of  the  values  of  the 
eight  independent  parameters  covered  by  the  correspondingly  numbered 
table  found  at  the  end  of  this  Memorandum.  If  the  key  indicates  (by 
an  x  in  the  parameter  matrix)  that  a  solution  is  available,  then  the 
numerical  values  are  given  in  the  correspondingly  numbered  table.  The 
key  also  provide:  the  reader  with  a  graphic  display  of  the  solutions 
available  in  the  neighborhood  ;f  his  range  of  interest.  Inner  solu¬ 
tions  for  the  limit  b  -*  »  are  given  in  Table  7  which  is  too  short  to 
warrant  a  separate  key. 

Approximately  a  sixth  of  the  solutions  listed  have  been  reported 
by  other  authors.  All  values  listed  have  been  computed,  or  recomputed, 
according  to  the  numerical  procedures  described  in  Section  IV.  Some 
differences  exist  between  the  present  values  ai  those  of  earlier  in¬ 
vestigators,  but  the  present  values  are  believed  to  be  accurate  with¬ 
in  the  limits  prescribed  in  the  succeeding  three  sections  (generally, 
correct  to  four  significant  figures). 

The  most  extensive  of  the  earlier  tabulations  of  solutions  to  the 


similarity  equations  foi  a  laminar  boundary  layer  of  a  perfect  gas  may 

be  found  in  the  works  of  Beckwith, Beckwith  and  Cohen, ^  Li  and 

Nagaraatsu,  Cohen  and  Reshotko,  Reshotko  and  Beckwith, and 

Dewey. Solutions  of  high  accuracy  have  been  obtained  by  Smith^^*^^ 

and  others.  Early  contributions  to  the  understanding  of  the  role  of 

( 131 

fluid  properties  may  be  found  in  the  works  of  Busemann,  ' 


Karma n  and 
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Tsien,^^  Crocco,^^  Young  and  Janssen,  and  var.  Driest. 

Emmons  and  Leigh''  ’  report  a  large  number  of  solutions  with  the  sur¬ 
face  mass-transfer  parameter  f^  other  than  zero.  Several  recent  com- 
pendiums^^  t0  on  boundary- layer  theory  may  be  consulted  for  a  more 
complete  description  of  previous  contributions. 

Since  this  Memorandum  is  restricted  to  positive  values  of  the 

pressure-gradient  parameter  3  and  values  f  £  0  corresponding  to  mass 

w 

injection  into  the  boundary  layer,  a  list  of  sources  is  offered  in 

( 18) 

which  solutions  for  3  <  0  and  f  >  0  may  be  found.  Emmons  and  Leighv  ' 

w 

have  obtained  solutions  for  3=0  and  values  of  f  /T  equal  to  10,  6,  , 

6,  3,  2.1,  2,  1.5,  1.4,  1*3,  ...  .5,  >45,  .40,  ...  0.1,  and  0.05,  as¬ 
suming  that  Pr  -  1  and  the  vi.'.cosity-temperature  law  is  linear.  Suc¬ 
tion  results  are  reported  also  by  Spalding  and  Evans, Pretsch,^^ 
Watson, Eckert,  Donoughe,  and  Moore,  Thwaites,^^  Schlichting 

and  Bussman,^^  Mangier, Schaefer, Enmons  and  Leigh, and 
(31) 

Koh  and  Hartnett.  ' 

Negative  values  of  the  pressure-gradient  parameter  3  have  been 
considered  by  Hartree,^^  Stewartson,  (^1*33)  Smith, Cohen  and 
Reshotko,^)  Beckwith, and  Hufen  and  Wuest.^^  Solutions  fcr  nega  ¬ 
tive  values  of  3  differ  from  their  counterparts  for  positive  3  in  two 
ways.  First,  two  sets  of  "proper"  solutions  exist  for  each  value  of 
3;  second,  every  value  of  the  wall  derivative  f^  and  its  corresponding 
value  of  0j,  etc.,  will  satisfy  the  boundary  conditions  f  '  =  0  at  T|  =  0 
and  f '  —  1  as  T)  -*  <»,  The  "proper"  solution  must,  therefore,  be  defined 
as  that  solution  for  whicu  t'  approaches  unity  most  rapidly  from  above 
(see  Refs.  3,  8,  and  21). 
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Key  1 

Sirailai  Solutions  for  u>  ■  Pr  ■  1 
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Key  1 ,  cont . 
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Key  2 


Similar  Solutions  for  a  Power-Lav:  Viscosity  Relation 

Pr  -  0.7,  f  -  0 
w 


Key  2,  cont 
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Similar  Solutions  for  Pr  =  0. 7,  t  *  1,  f  4  0 
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Key 


5 


(a) 


Similar  Solutions  for  a  Sutherland  Viscosity — Temperature  Relation, 


Pr  e  0*7,  t^  c  1.0 


a 

Two  sets  of  answers  are  given  for  each  case.  The  first  set 
represents  the  Sutherland  solution  for  the  indicated  value  of  s,  while 
the  second  set  represents  the  solution  using  a  power-law  viscosity- 
temperature  relation  with  u  e  ^  as  defined  by  Eq .  (53b)* 


Solution  independent  of  s. 
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Key  6 

Outer  Solutions  for  f3  -*  ® 
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IV.  NUMERICAL  INTEGRATION  PROCEDURE 

The  program  for  the  numerical  integration  of  the  system  of  equa¬ 
tions  (Eqs.  (40)  to  (43))  was  written  at  RAND  in  FORTRAN  for  the  IBM 
7044.  The  system  was  treated  as  a  two-point  boundary-value  problem, 
and  the  Runge-Kutta  method  was  employed  for  the  numerical  integration. 
The  sequence  of  operations  is  shown  schematically  in  Fig.  2.  Inputs 
to  the  program  are  values  for  the  eight  parameters  and  initial  guesses 
for  the  four  wall  derivatives  f"(0),  G^O),  g^O),  and  z'(0).  The  so¬ 
lutions  are  then  separated  into  two  categories:  ordinary  boundary- 
layer  solutions  and  adiabatic  wall  solutions.  The  latter  requires  a 
subprogram  that  searches  for  and  obtains  solutions  in  which  the  adia¬ 
batic  wall  condition  is  met.  In  either  category,  the  Runge-Kutta  meth 
od  is  used  to  integrate  the  equation  to  a  given  value  of  ’lTOax*  0n  the 
basis  of  experience  ...  is  value  was  selected  to  yield  an  acceptable 
asymptotic  solution  within  the  limits  of  the  variable  values  at  1]^  . 

The  integration  stepsize  was  varied  to  insure  the  proper  accuracy,  and 

the  values  of  the  variables  at  Tl  were  held  to  within  10  , 

max 

When  the  integration  procedure  is  completed,  the  final  values  at 
“0  are  coraPare^  with  the  required  values  and  the  Newton-Raphson 
scheme  is  used  to  recalculate  the  initial  condition.  The  procedure  is 
as  follows: 


Let 


INPUT  DATA: 
PARAMETERS 
&  INITIAL 
CONDITION 
GUESSES 


COMPUTATION 
OF  TW  BASED 
ON  INPUTS 


TAW  CASE 


COMPUTE 

CONSTANT 

EXPRESSIONS 


RUNGE-KUTTA 

INTEGRATION 

INTEGRATE 


/IS\.  YES  YES  YES  ^^IS^^ 

9  b max)  =  1  fomax)  =  1  ?>— <G (»max)  =  1  ?>— <['(w)  =  1  ?, 


TAW  CASE 


RECOMPUTE 

TW 


IS  NEW  TW 
sf  OLD  TW'^ 


ERROR- 

CORRECTION 

ROUTINE 


PUNCH  CASE 
DATA  &  PRINT 
PROFILE 


IF  THE  INITIAL  GUESSES 
ARE  SUCH  THAT  THE 
INTEGRATION  CANNOT 
BE  COMPLETED  THE 
PROGRAM  ALTERS  F"  (0) 
&  BEGINS  AGAIN 


Fig.  2  --  Pro.  ram  flow  chart. 
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t  =  e[f"CO),  g'(o),  z'(o),  e '(0) ]  =  f'OT^)  -  i  (82) 

3  -  3[f"(0),  g  (0) ,  z'(0),  e  '(0)1  =  g(\ax)  '  1  (83) 

*>  =  $[f"(0),  g  '(0),  z  /(0) ,  e  '(0)]  =  Q^\ax)  -  1  (84) 

M  =  w[f"(0),  g  ;(0) ,  z'(0),  e'(0)1  =  zOT^)  -  1  (85) 


['<U» 
8  IU 

6  <w 

J=  <WJ 


and  S  = 


r  fw(o) 


g'(0) 

e'(o) 

z\0) 


so  that  -e  «  -e(S)  =  B(S)  -  1  =  AB.  It  is  required  that  *  =>  *(£,  3,,8,  H)  ->  0. 
Now  d-*  =  dB  =  M  •  dS,  where 


dz  '(0) 

5gb 
dz  '  (0) 

56b 

5z  '(0) 


r  sfb 

< 

Sfb 

of  (0) 

dg  '(0) 

56 '(o) 

I  5gb 

3gb 

dg 

6b 

1  „ 

j  5f  (0) 

Sg'(0) 

56 '(0) 

S6b 

seb 

"b 

df"(0) 

ag'(o) 

56 '(0) 

Szb 

dzb 

9zb 

5f"(0) 

dg  '(0) 

56 '(0) 

-V  I 

'(o)_ 


The  integration  is  performed  using  the  initial  guesses  for  S. 
This  gives  B  and  «q(S).  The  initial  guesses  are  then  perturbed  as 
follows : 
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£[f  '(0)  +  Af  ^(0) ,  g'(0),  9'(0),  z'(0)]  (88) 

3[f"(0),  g'(0)  +  Ag'(0).  9  /(0) ,  z'(0)]  (89) 

*[f"(0),  g  ;(0) ,  9'(0)  +  A9'(0),  z'(0)]  (90) 

W[f"(0),  g  /(0)  ,  9'(0),  z'(0)  +  Az  /(0)  ]  (91) 


These  integrations  give  the  columns  of  M.  The  new  guesses  for  the  ini¬ 
tial  condition  are  corrected  by  calculating: 

AS  =  -M-1  -e  (92) 

o 


It  was  established  that  convergence  was  very  difficult  for  cases  involv¬ 
ing  high  &  (p  ^  2).  In  these  cases,  the  10  ^  limit  accuracy  was  relaxed, 
and  this  is  indicated  in  each  table  as  needed. 


1 
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V.  APPLICATIONS  OF  THE  CONCEPT  OF  LOCAL  SIMILARITY 

GENERAL  DISCUSSION 

The  numerical  results  tabulated  in  this  Memorandum  are  rigorously 
applicable  only  when  the  numerous  similarity  requirements  listed  in 
Section  II  are  satisfied.  In  considering  the  diverse  applications  of 
compressible  laminar  boundary-layer  theory,  it  is  a  rare  occurrence 
indeed  when  all  of  these  conditions  are  met.  A  most  important  ques¬ 
tion  then  arises:  What  approach  should  be  used  in  predicting  the  be¬ 
havior  of  the  nonsimilar  laminar  boundary  layer? 

The  copious  literature  relating  to  this  question  offers  four  ba¬ 
sic  types  of  approach.  The  first  is  to  abandon  the  similar  solutions 
entirely  and  adopt  approximate  techniques  such  as  integral  and  series 

solutions  containing  free  parameters.  The  most  successful  of  the  in- 

(35) 

tegral  approaches  appears  to  be  that  developed  by  Tani.  The  trans- 

(36) 

cendental  approximation  proposed  by  Hanson  and  Richardson'  '  and  the 

(37) 

"improved  approxiiu  cion"  technique  of  Yang'  '  also  appear  very  pro¬ 
mising.  Related  to  the  integral  methods  is  the  powerful  "strip  meth- 

(38) 

od"  proposed  by  Pallone,  which  follows  closely  the  inviscid  flow 

(39) 

integral  method  of  Belotserkovski. '  ' 

A  second  type  of  nonsimilar  calculation  employs  a  strictly  numer¬ 
ical  approach.  The  complete  nonsirailar  boundary- layer  equations  are 
used  and  a  new  set  of  calculations  is  performed  for  each  particular 
problem.  Examples  of  this  approach  may  be  found  in  the  works  of  Smith 
and  Clutter and  Fl:igge-Lotz  and  Baxter.  Although  numeri¬ 

cally  satisfying,  these  calculations  are  extremely  expensive  and  in¬ 
tractable  to  generalization. 
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(43) 

The  third  type  of  approach,  suggested  by  Lees,  has  been  most 
successful  in  capturing  the  spirit  of  the  use  of  similar  solutions  in 
situations  where  exact  similarity  does  not  exist.  He  observed  that  un¬ 
der  certain  circumstances,  notably  when  there  is  a  highly-cooled  body 
in  hypersonic  flow,  the  local  pressure-gradient  parameter,  3,  had  a 
negligible  effect  on  the  heat  transfer  to  the  surface.  Many  elabora¬ 
tions  of  this  approach  have  been  proposed  to  improve  Lees'  simple  re¬ 
sult  to  provide  more  accurate  numerical  estimates  of  heat  transfer, 

(22) 

skin  friction,  and  boundary- layer  thickness.  Moore  gives  a  lucid 
summary  of  one  group  of  these  results.  Additional  ideas  for  modifying 

local  similarity  not  discussed  by  Moore  may  be  found  in  the  papers  of 

(44)  (45)  (6) 

Smith,  '  Kemp,  Rose,  and  Detra,  and  Beckwith  and  Cohen.  ' 

This  third  approach  is  saddled  with  one  significant  difficulty:  it 
is  necessary  to  make  one  or  more  implicit  ad  hoc  approximations  regard¬ 
ing  the  contribution  of  the  nonsimilar  terms  in  the  complete  boundary- 
layer  equations.  In  each  of  the  papers  cited  in  the  previous  paragraph, 
the  question  is  not  “Are  similar  solutions  applicable?"  but  rather  “Which 
similar  solution  should  be  used?"  A  number  of  methods  have  been  proposed 
for  choosing  the  similar  solutions  most  appropriate  to  the  local  inviscid 
flow  conditions,  local  wall  tempprature,  and  boundary- layer  history.  In 
the  context  of  the  present  discussion,  this  question  necessitates  the  ju¬ 
dicious  choice  of  values  for  the  eight  similarity  parameters.  These  val¬ 
ues  are  usually  determined  (e.g.,  Ref.  6)  by  satisfying  one  or  more  inte¬ 
gral  conservation  equations  exactly  using  assumed  profiles  obtained  from 
similar  solutions.  Such  procedures  are  closely  related  to  the  integral 
method  of  Thwaites^^  and  to  the  more  recent  use  by  Lees  and  Reeves'"4^ 

(33) 

of  a  family  of  similarity  profiles  generated  by  Stewartson. 
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The  fourth,  and  in  many  ways  most  satisfying,  type  of  approach  to 

nonsimilar  boundary- layer  calculations  may  be  traced  to  the  work  of 
(19) 

Meksyn.  His  underlying  premise  is  very  powerful:  if  the  boundary 

layer  is  at  all  times  very  nearly  described  by  a  similar  solution,  then 
the  direct  effects  of  the  nonsimilar  terms  may  be  calculated  by  asymp¬ 
totically  expanding  the  full  boundary- layer  equations  in  terms  of  small 
parameters  which  measure  the  departure  of  the  solutions  from  similarity. 
In  this  way,  the  accuracy  of  local  similarity  methods  is  explicitly  de¬ 
termined  by  using  the  full  nonsimilar  equations.  We  shall  demonstrate 
shortly  that  the  linearized  equations  governing  the  departure  from  sim¬ 
ilarity  depend  only  on  the  local  similarity  parameters  and,  consequently, 
need  be  computed  only  once. 

ASYMPTOTIC  EXPANSION  OF  THE  BOUNDARY - LAYER  EQUATIONS 

For  purposes  of  illustration,  we  shall  limit  our  consideration  to 
the  "incompressible"  momentum  equation 

+  ££Htl  +  W5)ri  '  £^  *  2![Vs-  -  £!£t,T|1  <«> 

where  §,  T|,  and  3(?)  are  as  defined  in  Section  II  and  the  subscripts  T| 
and  |  denote  derivatives.  Equation  (93)  may  be  obtained  from  the  gen¬ 
eral  equations  by  assuming  Pr  =»  u)  =  tg  =  t  =  1.  We  also  assume  that 

f  =  0,  making  the  three  boundary  conditions  for  Eq.  (93) 
w 

f(?,o)  =  o 
f^G.o)  =  o 

f^(l,°°)  =  i 


(94) 
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(49) 

Merk  was  the  first  to  expand  the  complete  nonsimilar  momentum 
equation  in  terms  of  a  small  parameter.  Subsequently  Bush^"*^  pointed 
out  that  Merk's  derivation  neglected  important  terms  in  the  correction 
equation.  The  remainder  of  this  section  will  be  concerned  with  Bush's 
equations  and  their  approximate  solution. 

In  the  spirit  of  Meksyn's  approximation,  we  look  for  solutions  to 
Eq.  (93)  when  the  right-hand  side  is  small.  The  key  to  an  appropriate 
expansion  is  the  inversion  introduced  by  Merk.  We  change  variables 
from  [§,T|;  0(§)]  to  [0,T|j  §(&)]  so  that  the  streamwise  momentum  equa¬ 
tion  becomec 

fTyTF)  +  ffT)Tl  +  ^  “  fT^  "  e<e)C-Tjfpf7j  -  (95> 

£(0,0)  -  ^(0,0)  -  0  ;  ^(0,®)  -  1  (96) 

where 

e(0)  -  2§e'(S)  -  25(p)/5#(0)  (97) 

If  e  is  zero,  the  momentum  equation  reduces  to  the  Falkner-Skan  simi¬ 
larity  equation  with  0  as  a  single  parameter.  For  small  e  we  may  per¬ 
form  an  asymptotic  expansion  of  f(3,T|)  of  the  form 

£(0,11)  -  fo(0,T»  +  s(p)f1(p,Tl)  +  ....  (98) 

Then  the  derivatives  f_  and  fQ  are 

M  P 

£„  ■  <£0>n  +  e(BH£i>n  + 


0  0  0 


(99) 
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f8  *  +  +  ••• 

+  e  /(g)f1  +  ...  (ICO) 

As  Bush  pointed  out,  the  term  e  /(P)  is,  in  general,  of  order  unity  and 
may  be  expressed  as 

1  d 

e'(P)  - - —  [25&/(§)]  =  2[1  +  E(8)]  (101) 

f3'(§)  dS 

where 

E(p)  =  SP"(S)/P'U)  =  -§(P)  §"0)/C§'0)l2  (102) 

Substitution  of  the  asymptotic  sequence  [Eq.  (98)]  into  the  momen¬ 
tum  equation  and  boundary  conditions  produces  a  hierarchy  of  equations, 

the  first  two  of  which  are  (primes  on  f  and  x,  denote  differentiation 

o  1 

with  respect  to  71)  : 

Order  Unity 


f"  +  f  f"  +  8[1  -  (f  ')2J 

0  O  O  0 


Order  e 


+  fofl  -  Alfofl  +  A2fofl  =  *<M) 

^(M)  *  f'(P,0)  -  f^P,-)  *  o 


fo(M)  =■  fJO,0)  =  0  ;  fQ'(P,»)  =  1 


(103) 


(104) 
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In  Eqs.  (104)  the  terms  (A^  A  ,  $)  are 

Aj^  -  2  +  28  +  2E 

A2  »  3  +  2E 

»(M)  -  *,'<*„'>(!  -  <Ve£: 

Note  that  two  independent  parameters  (8,E)  appear  in  the  first-order 
equation. 

Equations  (103)  represent  t\z  similar  solutions  of  Falkner  and 
Skan.  Equations  (104)  represent  the  first  correction  f^  to  the  ve¬ 

locity  profile  which  arises  from  nonsimilar  terms.  For  example,  the 
skin-friction  derivative  f"(8,0)  is  expressed  as 

^'(0,0)  -  f£(M)  +  «lj(0,u)  +  ...  (106) 

where  fQ(0>O)  is  the  local  similarity  solution  corresponding  to  the  lo¬ 
cal  value  of  0,  and  f"(P,0)  is  the  correction  obtained  as  a  solution  of 
Eqs.  (104).  It  is  apparent  that  the  correction  will  be  of  order  e  as 
long  as  1^(8, 0)  ic  of  the  same  order  ap  f"(8,0)  and  e  is  smal1. 

DPTERMI NATI ON  OF  f^lft.Tp  BY  SUCCESSIVE  APPROXIMATIONS 

We  proceed  to  a  consideration  of  the  first-order  equations  for 
^(M),  Eqs.  (104).  The  differential  equation  is  linear  with  homoge¬ 
neous  boundary  conditions  and  may  be  solved  numerically.  The  primary 
difficulty  arises  in  computing  the  inhomogeneous  term  $(8,1))  which  con¬ 
tains  derivatives  of  fQ  with  respect  to  both  7[  and  0.  This  is  a  diffi¬ 
cult  term  to  calculate  numerically  because  a  number  of  similarity 
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solutions  in  the  neighborhood  of  8  must  be  known  with  high  precision. 

The  exact  numerical  calculation  of  f^  appears  possible  but  has  not  been 
attempted. 

The  technique  adopted  here  is  to  substitute  a  transcendental  ap¬ 
proximation  for  fo(M)  which  allows  the  coefficients  of  the  differen¬ 
tial  •,uation  and  the  forcing  function  $0,11)  to  be  expressed  in  terms 

(52) 

of  known  functions.  Following  an  earlier  paper  by  Bush,  we  repre¬ 
sent  f  by  the  relation 
o 

f0'O,Tl)  =  erf(t)  (107) 

where 

t  -  aT|  (108) 

aO>  -  -  j 1  +  e 

This  relation  is  found  to  be  in  excellent  agreement  with  exact  solutions 
for  fQ  and  appears  quite  adequate  for  our  present  purposes.  It  is  con¬ 
venient  to  transform  coordinates  from  0,11)  to  (a,t)  and  define  a  new 
dependent  variable  g(a,t)  according  to  the  relation^ 


3 

3 

3 

/  t  3  3  \ 

■  ■  ■  s 

a  — 

, 

*  a'O) - +  — 

(110) 

311 

3t 

O 

'  a  3t  3a  ' 

f^M) 

=  X(a)g(a,t) 

(111) 

2  /  da  \  0. 23084 

X(a) 

=  ~2~Z\  — )  -  3 

(112) 

a  Vtt  '  d&  '  a 
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Substitution  of  those  transformations  into  Eq.  (104)  results  in  a 
linear  differential  equation  for  g(a,t)  which  contains  derivatives  of 
g  with  respect  to  t  up  to  third  order.  The  boundary  conditions  for 
this  equation  are  (primes  denote  different! stion  with  respect  to  t) : 


g(a, 0)  =*  g  /(a,0)  =  g'(a,»)  =  0  (113) 

The  differential  equation  for  g(a,t)  is  now  integrated  formally 
three  times  with  respect  to  t,  using  the  boundary  conditions  given  by 
Eq.  (113).  The  resulting  integral  equation  is 


g(a, t)  *  —  g  (a,0)  -  J  J1(t)g(a,t)  dt  +  A3  [  JjQ(t)g(a,t)  dt 


-A4  J  J  f  J_1(t)g(a,t)  dt  +  t1(t) 


(114) 


The  coefficients  A^(a)  and  A^(a)  are  given  by 

A3(a)  =  4  +  2E  +  23  ;  A4(a)  =■  6  +  4E  +  23  (115) 

and  the  terms  Jn(t)  are  the  n^  integrals  of  the  error  function: 


2  2 

n  -  -UP  J.x(x)  -  y=  e‘x 


J  (x)  ■  erf(x) 

o 


:  •  Jirfl(x) 


1  +  n 


xj  (x)  +  -  J  .  (x) 
n'  2  n-1 


x 

n  .'VtT 


(116) 


The  term  ^(t)  is  defined  by 
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i  i  i 

*,(t)  =  —  -  j,(V5t)  -  -  j,(t)  -  -  jQ(t)j .(t)  +  -  *_(t)  (117) 

1  2V2  2  2  2  8  °  1  4  2 


where 


*2(t) 


1  1 

-  j  <*K>  -  -  J0(t> 
V2  2 


+ 


Vn 

—  J  ,(c) 

2 


1 


(118) 


A  method  of  successive  approximations  is  now  applied  to  Eq.  (114), 

V 

substituting  trial  functions  g(a,t)  for  g(a,t)  in  the  three  integrals 
and  continuing  until  a  trial  function  g  is  found  that  agrees  satisfac¬ 
torily  with  the  function  g  computed  from  the  integral  equation.  This 
scheme  differs  from  Picard's  method  in  that  the  sequence  of  trial 
functions  g  used  in  the  integrals  are  suitably-chosen  integrable  func¬ 
tions  of  t  rather  than  the  functions  g(a,t)  obtained  in  the  previous 
iteration  step.  Picard's  method  converges  absolutely,  but  in  practice 
it  usually  cannot  be  continued  analytically  beyond  one  or  two  itera¬ 
tions.  In  the  present  scheme,  convergence  depends  on  the  choice  of 
trial  functions  g  but  the  integrals  may  be  evaluated  in  closed  form. 

Evaluation  of  the  Velocicy  Profile 

In  this  analysis  based  on  the  simplified  momentum  equation  [Eq. 
(93)]  we  are  most  interested  in  the  correction  f''(8,0)  to  the  velocity 
gradient  at  the  wall.  It  is  therefore  more  convenient  to  work  with 
the  first  derivative  of  Eq.  (114),  which  is  (after  some  rearrangement) 


g'(a,t) 


tg"(a,0)  -  A^J j(t)g(a,t)  +  Afe 


dt 


+  A4  J  {  J0(t)g'(a,t)  dt  +  *3(t) 


(119) 
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vliere 

A.  *  3  +  2E  ;  A,  -  2  +  2E  (120) 

J  o 

1  1  Vn 

*,(t)  -  -  J.(V2t)  -  -  J.(t)  -  —  j‘(t)  (121) 

3  2  1  2  1  8  ° 

In  following  the  technique  of  successive  approximations,  the  terms 
g  and  g'  on  the  right-hand  side  of  Eq.  (119)  must  be  replaced  by  the 
trial  function  g  and  its  derivative  g#.  The  behavior  of  g  and  g '  i  ay 
be  inferred  from  the  boundary  conditions  g(a,0)  ■  g'f.a.O)  =  g,(a,or)  ■  0 
and  is  sketched  in  Fig.  3« 


t 


Fig.  3  --  Behavior  of  g  and  g'. 
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The  behavior  of  f^(t)  and  the  two  integrals  appearing  in  Eq.  (119)  may 
be  inferred  from  the  behavior  of  g'  and  the  quantities  Jn(t). 

Suppose  we  decompose  the  trial  function  g(a,t)  into  two  parts,  so 

that 

g(a,t)  -  C(a)H(t)  (122) 

We  then  define  n(t)  to  be  unity  as  t  -»  •,  making  g(a,“)  *■  C(a)  4  0. 

The  boundary  conditions  g(a,0)  ■  0  and  g^a.O)  ■>  0  are  automatically 
satisfied  for  all  suitable  trial  functions  g(a,t);  tue  boundary  condi¬ 
tion  g'(a,«)  =*  0  serves  to  evaluate  both  C(a)  and  g'(a,0).  The  terms 
appearing  on  the  right-hand  side  of  Eq.  (119)  have  the  following  asymp¬ 
totic  forms: 


lim 

t-*® 


[j 1(t)g(a,t) ] 


lim 

t-<® 


J1(t)g/(a,t)  dtj 


c(a>V1 


lim 

t-® 


[U 


JQ(t)g/(a,t)  dtj 


C(a)  [y2t  -  v3] 


lim 


1  _  i/n 

-  (V2  -  l)t  - 

2  8 


(123) 

(124) 

(125) 

(126) 


The  terms  (Y^,  Y2>  ^  )  are  numerical  conjtants  which  art  determined  by 
the  choice  of  the  trial  function  n(t). 

Substituting  Eqs.  (123)  to  (126)  into  Eq.  (1.19)  and  applying  the 
boundary  condition  g^a,®)  a  0,  we  obtain  the  following  formulas  for 
C(a)  and  g"(0): 
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g /(a,0)  =  C(a) [(1  +  As)  -  A4Y2]  -  -  (VI  -  1) 


(127) 


(128) 


We  have  chosen  several  functions  n(t.)  and  used  Eqs.  (127)  and 
(128)  to  determine  C(a)  and  g"(a,0).  The  integrals  appearing  in  Eq. 
(119)  have  been  evaluated  in  closed  form  and  the  profiles  g'(a,t)  have 
been  determined  for  each  trial  function. 

Products  and  sums  of  the  functions  J^(t)  are  useful  in  generating 
successively  mere  accurate  trial  functions  n(t)  because  the  integrals 
are  easily  evaluated.  Typical  examples  of  n(t)  are: 


K(t) 


YTT  | 

--J_l(t)| 


(VI  -  1)  2 


(t) 


(129) 


and  there  are  others. 

Weighted  sums  of  the  functions  listed  above  were  also  used  to  ef¬ 
fect  a  close  fit  to  g/(a,t).  The  computed  profiles  g'(a,t)  differed 
in  magnitude  with  the  different  functions  a(t),  but  both  the  qualita¬ 
tive  behavior  and  the  quantity  g^ajO)  were  relatively  insensitive  to 
the  form  of  n(t). 


-55- 


For  P  =  E  *  0,  g'^a.t)  has  a  maximum  near  t  =  0.84  and  decreases 
to  about  12  percent  of  its  maximum  value  at  t  =  2.0.  The  derivative 
g  (a,0)  is  approximately  -0.055.  For  3=0,  the  term  a(3)  computed 
from  2q.  (109)  is  1/2,  and 

3  °  E  -  0 

fi<.0,0)  =  a2x(a)g"(a,0) 


-.025 


The  Blasius  solution  to  Eqs.  (103)  for  3  -“  0  is  f"(0,0)  =  0.46960,  and 
the  final  estimate  of  the  correction  to  the  s.cin  fraction  coefficient 
for  3  *  E  =  0  is 

f"(0,0)  a:  f"(0,l)[l  -  0.053  e  +  ...  ]  (130) 

The  physical  interpretation  of  this  result  agrees  qualitatively 
with  the  sign  of  the  correction-  If  e  >  0,  3  is  increasing  and  the 
local  similarity  value,  f^(3,0),  would  neglect  "relaxation"  effects 
and  overestimate  the  shear  at  the  wall.  The  correction  reduces  the 
skin-friction  value  by  an  amount  proportional  to  e.  The  small  magni¬ 
tude  (.053)  of  the  correction  term  is  somewhat  surprising,  but  this 
value  is  probably  accurate  to  within  about  ±  20  percent.  Bush  esti¬ 
mated  the  correction  terr  to  be  0.204  but  his  result  was  obt'ined  by 
asymptotically  expanding  an  approximate  solution  rather  than  obtaining 
an  ap£ coximate  solution  to  the  exact  first-order  equation. 


* 
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VI.  SOLUTIONS  FOR  LARGE  VALUES  OF  THE  PRES SURE -GRADIENT  PARAMETER  B 


In  using  the  coueept  of  local  similarity  in  highly  accelerated 
flows,  it  is  useful  to  have  solutions  of  the  laminar  boundary -.layer 
equations  for  large  values  of  P.  The  mathematical  difficultieo  en¬ 
countered  in  the  li  it  P  »  1  may  be  illustrated  by  rewriting  the 
streauiwise  momentum  equatin .  [Eq.  (40)]  in  the  form 


^  Of")  7  +  ff"]  -  |f /2  -  —  [(1  -  tw)0  -  (1  -  ts)g2  +  tjj  -  0  (131) 


Here  X  13  redefined  as  the  density-viscosity  racio  at  the  wall  (p^/p^p.^) 


In  the  limit  P  -•  »,  Eq.  (131)  reduces  from  third  to  first  order: 


lim 

P-“ 


(  1  f 

'  '  17  l(1  • 


)e  -  (i  -  ts)g 


2  ♦  vf2 


(132) 


The  transverse-momentum  equation  and  the  energy  equation  remain  of  sec- 

(53) 

ond  order.  As  originally  pointed  out  by  Coles  and  as  later  elabor¬ 
ated  upon  by  Beckwith  and  Cohen., this  leads  to  a  singular  perturba¬ 
tion  problem  in  which  the  thickness  of  the  velocity  layer  is  of  order 
-1  2 

P  with  respect  to  the  total  enthalpy  layer  and  transverse  ve.oci  ty 
layers  of  order  unity.  The  method  of  solution  is  similar  to  that  em¬ 
ployed  in  deriving  a  uniformly  valid  approximation  to  the  Navier-Stokea 

r  f55) 

equations  in  the  limit  of  large  Reynolds  number  [see  Kaplun, '  Lager- 
strcms  and  Cole,^^  and  the  recent  book  by  Van  Dyke^2  ]. 

*Discussicn  of  this  problem  also  appears  in  an  abbreviated  version 
(54) 

in  Lagerstrom's  article. 


* 
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THE  OUTER  LIMIT  EQUATIONS 

Since  the  mathematics justification  of  this  singular  perturbation 

solution  (more  popularly  called  an  inner-  and  outer-expansion  procedure) 

(5^) 

has  been  discussed  in  detail  by  C  ies,v  '  cur  purposes  will  be  served 
by  a  cursory  development  of  the  governing  equations.  The  present 
analysis  extends  the  work  of  Beckwith  and  Cohen^  to  include  a  power- 
law  temperature-viscosity  relation,  a  constant  but  nonunit  Prandtl  num¬ 
ber,  and  the  mass  transfer  at  the  surface. 

Assume  that  appropriate  outer  representations  of  the  dependent 
variables  f,  8,  g  are  of  the  following  forms: 


1 

f  =  f  +  —  f ..  +  ... 

»  v?  1 


0  =  8+  —  6.  +  ... 
°  Vb  1 


8  -  80  +  -  »1  +  ~. 


(133) 


Substituting  these  representations  into  the  momentum  equations  and  the 

-1/2  * 

energy  equations  and  dropping  all  terms  of  order  0  and  smaller, 
the  following  "outer"  equations  are  obtained: 


The  outer  limit  equations  are  properly  obtained  by  applying  rhe 
limit  0  -»  oo  to  the  full  equations  expressed  in  outer  variables,  with 
7]  held  fixed.  This  gives  identical  results  to  those  cited  here. 
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<<>  /  +  f0*o 


(136) 


The  appropriate  boundary  conditions  are  found  (a)  by  requiring 
that  the  outer  equations  satisfy  the  exact  ou"er  boundary  conditions, 
and  (t)  by  exact  matching  r *  the  inner  representations  with  the  outer. 
The  results  may  be  written 

f  (0)  -  f  ;  6  (0)  -  g  (0)  =  0  ;  0  (»)  -  g  (»)  *  1  (137) 

O  W  o  O  O  O 


Only  one  boundary  condition  on  f  may  be  satisfied  by  the  outer 

equations,  because  the  outer  limit  equation  for  fQ  reduces  to  first 

order.  The  exact  boundary  condition  f  ;  -*  1  as  1)  -»  <*>  is  automatically 

satisfied  by  Eq.  (134).  The  no-slip  condition,  f ;  -*  0  as  T|  -*  0,  must 

be  satisfied  by  an  inner  solution  which  is  valid  in  a  region  of  extent 

“  i  /2 

0  with  respect  to  the  scale  of  the  outer  solution. 

In  this  approximation,  the  density-viscosity  ratio  X  =  (pn/p  P-  ) 

w  w 


in  the  outer  layer  becomes 
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X 


(138) 


where 


g2 


2 

cos 


A  +  sin 


2 


A 


As  T|  -•  0,  X  —  [(1  -  a.)  ft  l^-1  because  f '(0)  =*  (t  /t  )^2  and  f  '(0)  is 

^  8  O  W  S  O 

* 

nonzero  in  general-  The  energy  and  transverse  motnei  ;.u<n  equations  are 
coupled  through  the  implicit  appearance  of  Doth  9q  and  gQ  in  Eqs.  (135) 
a  d  (136).  If  Pr  *»  1,  0q  -  gQ,  and  the  number  of  coupled  ordinary  dif¬ 
ferential  equations  is  reduced  from  three  to  two.  Setting  both  the  vis¬ 
cosity-temperature  exponent  and  the  Prandtl  number  Pr  equal  to  unity 
reproduces  the  equations  of  Eeckwith  and  Cohen.^  If  =  0  (i.e. , 

the  local  Mach  number  is  zero),  t  =»  1,  and  the  transverse  momentum  and 

s 

energy  equations  are  again  uncoupled. 

One  interesting  case  was  pointed  out  by  Coles  and  we  extend  his  re¬ 
sult  to  generalized  compressible  flow.  Let  t  “  t  ■  1,  A  ■  0  so  that 

w  s 

ct„  ■  (U2/2H  )(u  /u  )2  =  o-  •  Then  X  ™  (1  -  c1)a’"1  and  f '  1  for  all  T). 

/  oo  0  e  oo  1  i  o 

2  *1 

In  taking  this  limit  with  Pr  ^  1,  the  product  of  (fj  )  and  (1  -  t^) 
approaches  0j  so  that  the  energy  equition  becomes 


0"  +  0'cn  +  f  )  Pr  (1  -  aj^l  -  a, (1  -  Pr)]"1  =  0  (139) 

O  O  W  i  1 


It  should  be  noted  that  the  inner  and  outer  expansion  procedure 

breaks  down  in  the  limit  t  -*  0,  because  the  outer  solution  for  f  sat- 

w  o 

isfies  the  exact  boundary  conditions  for  the  complete  equations  and  the 
inner  solution  for  t’  is  simply  zero. 
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wt th  the  boundary  conditions 


T1  -  0  ,  0  (0)  -  0 

o 

T)  — »  00  ,  0^(oo)s»l 


(140) 


If  we  define  the  new  variable  x  and  the  constants  xo  and  W  by 

x  =  V7 ' 

f 
w 

x°  = 

w  =  [1  -  ax(  1  -  Pr)][Pr(l  -  ) 

then  Eqs.  (139)  and  (140)  are  satisfied  by  the  solution 

0  *  [erf(X/l^)  +  erf(x0/V2)]/[l  +  erf(xo/V2)]  ;  (142) 

For  xc  <  0,  we  note  the  identity 

erf(-X)  *  -erf(x)  (143) 

Equations  (134)  to  (136)  with  the  boundary  conditions  of  Eq.  (137) 
represent  a  two-point  boundary-value  problem  for  the  complete  solution 
of  the  outer  limit  equations.  The  derivatives  8^(0)  and  g '(0)  along 
with  the  integrals  1^  1^1),  1^2),  and  1^3)  evaluated  using  f  0  , 
and  gQ  in  place  of  f,  8,  and  g  are  given  in  Table  6.  The  integral  I 
is  identically  equal  to  zero.  Beckwith  and  Cohen<“6)  calculated  several 
of  these  quantities  for  the  case  of  Pr  -  u>  -  1,  f  .  o,  and  nonunit 


values  of  t  . 

s 


THE  INNER  LIMIT  EQUATIONS 


Inasmuch  as  the  order  o£  the  energy  and  transverse  momentum  equa¬ 
tions  is  not  reduced  in  taking  the  limit  g  -*  ®,  the  outer  equations 
[Eqs.  (134)  to  (136)]  represent  complete  solutions  for  the  total  en¬ 
thalpy  and  transverse  velocity  profiles  for  large  0.  The  terns  fQ  and 
2 

(fr  )',  which  appear  in  the  outer  equations,  differ  from  the  exact  so- 

2  - 1/2 

lutions  f  and  (£ '  )'  only  in  a  region  which  is  0  smaller  in  extent 

than  the  region  of  applicability  of  the  outer  equations.  Therefore, 

the  outer  equations  asymptotically  represent  the  complete  solutions 

for  9  l  id  g  as  g  -•  «,  and  in  this  limit,  the  inner  solutions  for  9  end 

g  are  identically  zero. 

The  no-slip  condition  f/(0)  *  0  is  satisfied  by  the  inner  limit 
equation  for  £;  the  inner  equations  for  9  and  g,  as  noted  previously, 
reduce  to  9  ■  g  ■  0.  To  examine  the  inner  streamwise-momentum  equa¬ 
tion,  it  is  necessary  to  introduce  the  new  independent  variable  and 
an  inner  representation  for  f: 

1 

T)  -  /Fl  J  f  -  f  +  ~  f_(f[)  +  ...  (144) 

Vg 

so  that 

f'(T»  "  f'CH)  -  —  (145) 

u 

e 

After  rearrangement,  the  inner  equation  for  f^  becomes 
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-1/2 

and  dropping  terms  of  order  &  and  smaller,  the  result  is 


<*V  ■  <fo> 


f ')2  +  — 
t 


(147) 


where 


I  a.  cos  A  „  »lm-l 

1  •  -h— 

w  J 


(148) 


The  boundary  conditions  on  are  found  from  the  exact  boundary  condi- 

ic 

tions  at  the  wall  and  the  matching  condition  that  the  inner  and  outer 
representations  of  f  agree  in  the  limit  V8  -*  00  with  large  but  fixed.’ 
The  following  boundary  conditions  for  the  inner  equation  result: 


fo'(0)  -  0  ; 


-  /yr 


(149) 


Since  Eqs.  (147)  and  (148)  involve  only  fj  and  not  tQ,  they  repre¬ 
sent  a  well  posed,  second-order,  two-point  boundary-value  problem.  A 
more  syumetrical  form  may  be  obtained  by  one  final  transformation: 


80<t) 


-v? 


t  -  (t  /t  ) 
W  8 


1/4 


(150) 


Then  the  problem  may  be  written 


(x.p #  +  i 


-  0 


(151) 


*The  matching  condition  applied  hare  is  elaborated  by  Van  Dyke: 
Inner  representation  of  "outer  representation"  -  Outer  representation 
of  "inner  representation". 
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with  the  boundary  conditions 

sQ(0)  -  0  ;  so(«)  «  1  (152) 


and  the  definitions 


X 


[1 


2,u)-l 
as  J 
o 


(153) 


a.  cos A  A/t 
1  s 


(154) 


Values  of  8^(0)  satisfying  Eqs.  (151)  to  (153)  are  given  in  Table 

7  for  several  values  of  tu,  a,  and  A.  For  the  special  case  of  a-  =  1, 

(53) 

Coles  pointed  out  that  an  analytic  solution  may  be  found  in  the  form 


80  =  1-3  sech2  (t/vT  +  tanh"1  V2/3) 


055) 


Using  either  the  analytic  solution  for  w  =  1  or  the  numerical  solutions 
for  w  j  1,  the  surface  skin  friction  derivative  f"  is  then  found  by  re¬ 
versing  the  previous  transformations 


lim  .// 

&-wo  w 


056) 


In  general  3^(0)  depends  on  the  three  parameters  a  ,  uu,  and  t  ;  for 

“  8 

^  *  1,  the  value  of  s^(0)  is  4/3. 

* 

Our  thanks  go  to  J.  Aroesty  of  The  RAND  Corporation  for  suggest¬ 
ing  a  simple  numerical  quadrature  technique  for  solving  Eq.  (151). 
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VII.  DISCUSSION 


The  large  number  of  the  solutions  in  Tables  1  to  7  (p.  78  ff.) 
suggests  many  possibilities  for  numerical  correlations  and  comparisons 
with  approximate  results.  Although  an  exhaustive  discussion  of  these 
possibilities  goes  beyond  our  present  purposes,  we  present  selected 
examples  in  Figs.  4  to  14  (p.  72  ff.)  of  the  use  of  the  present  solu¬ 
tions  in  understanding  the  influence  of  the  similarity  parameters  on 
heat,  mass,  and  momentum  transfer. 

The  local  skin-friction  and  heat-transfer  coefficients  C,  and  C, 

n  i 

are  defined  by  Eqs.  (70)  and  (73).  Following  Ref.  59,  they  may  be 

placed  in  a  more  symmetric  form  by  using  the  Reynolds  number 

Re  =  P  u  x/V  (157) 

e,x  e  e  e 

and  the  dimensionless  streamwise  coordinate 


S(p  p  u  r^x)  1 
w  w  e  k 


(f  V»VkJ  d*)  l».Wi,l'1 


(158) 


With  these  definitions,  the  quantities  and  become 


cos  A 


h  pi  \  1  / 2  /p  u  \ 

H  — 

\P  pi  /  \p  U  / 


e,x  e  e 


0'(1  -  t  ) 
w  w 


Pr (t  -  t  ) 

aw  w 


(159) 


2  cos  A 


2  ,  9  t 

cos  A  +  (g ')  sin“  A 
w 


1/2 


e,x  e  e 


(160) 
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In  Eqs.  (159)  and  (160),  only  the  terms  appearing  in  the  square 
brackets  depend  upon  the  similarity  parameters;  the  terms  preceding 
the  square  brackets  represent  the  external  flow  conditions  and  wa] 1 
temperature.  The  Reynolds  analogy  function 

j  -  2Ch/Cf  (161) 

is  simply  the  ratio  of  the  square  brackets  appearing  in  and  and 

consequently  depends  only  upon  the  similarity  parameters. 

In  Fig.  4,  the  expected  increase  in  the  skin-friction  coefficient 

with  increasing  pressure-gradient  parameter  8  is  shown.  The  relative 

increase  in  muc^  greater  for  higher  wall  temperatures. 

This  occurs  because  the  boundary  layer  is  thicker  and  the  response  to 

increasing  p  is  proportionately  higher. 

Figure  5  sh'iws  similar  behavior  for  the  heat- transfer  coefficient. 

We  have  used  the  wall-gradient  parameter  9/[(l  -  t  )/(t  -  t  )]  in 

W  W  flW  w 

(4)  . 

forming  the  heat-transfer  coefficient.  A3  shown  previously,  6^ 

•ff 

varies  greatly  with  P  for  t  >0.6  because  t  varies  with  P  if  o  ^  0. 

w  aw 

The  percentage  increase  here  is  less  than  for  the  skin  friction  but  the 

(43) 

trends  with  P  and  t  are  the  same.  Lees  argued  chat  the  high  den- 
w 

sity  in  the  boundary  layer  near  the  wall  for  low  wall  temperatures  "in¬ 
sulates"  the  wall  from  pressure-gradient  effects  and  the  influence  of 

* 

It  is  particularly  important  to  determine  the  proper  adiabatic 
wall  temperature  ir  cases  with  mass  injection  because  t  is  greatly 
decreased  by  injection  (see  Table  3).  The  assumption  Pr  *  1  predicts 
t  *1  under  all  conditions,  which  is  seriously  in  error  for  large 
injection  and  high  local  Mach  numbers. 
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0  on  is  greatly  reduced.  This  behavior  is  shown  with  decreasing 

values  of  t  . 

w 

The  next  two  figures  deal  with  the  Reynolds  analogy  function 

j  -  2C^/C^.  Li  and  Gr  showed  earlier  that  large  deviations 

from  unity  could  occur  in  a  hypersonic  boundary  layer  even  for 

t  «  1  and  p  ■  0. 
w 

Figure  6  demonstrates  that  j  decreases  with  increasing  P,  the  de¬ 
crease  being  faster  with  increasing  wall  temperature.  This  result  can 
be  easily  explained  by  examining  the  limiting  equations  for  P  -»  ®.  The 

heat-transfer  coefficient  C  approaches  an  asymptote  as  P  -»  whereas 

n 

1/2 

for  large  P  the  skin-friction  coefficient  increases  as  P  •  The  ap¬ 
proach  to  this  limitii^  behavior  is  more  rapid  with  larger  wall  temper¬ 
atures.  Similar  behavior  is  shown  in  Fig.  7  for  values  of  the  hyper¬ 
sonic  parameter  cj  -  0,  0.5,  and  1.0- 

Figure  8  shows  the  well  known  boundary- layer  property  that  heat 
transfer  is  decreased  by  mass  injection  at  the  surface.  The  heat- 
transfer  reduction  caused  by  a  given  value  of  f^  decreases  with  in¬ 
creasing  P.  The  reason  will  be  explained  below  in  the  discussion  on 
limiting  solutions  for  large  P. 

Finally,  the  sweep  angle  parameter,  t  ,  is  varied  in  Fig.  9.  When 

s 

the  sweep  angle  A  is  zero,  there  is  skin  friction  only  in  the  x-direction 
(see  Eq.  (160)).  As  the  sweep  angle  increases,  the  spanwise  term, 
g'(0)sin  A,  contributes  increasingly  to  the  skin-friction  coefficient. 

If  the  free-stream  direction  is  sufficiently  oblique,  the  x-component 
of  C^,  f"(0)  (ue/ue})  cos  A,  no  longer  exerts  an  appreciable  influence. 


The  quantity 
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{l  +  [g"(0)  sin  A/f"(0)  cos  a]^J 

plotted  in  Fig.  9  is  a  relative  measure  of  the  contribution  of  the  two 

skin-friction  components,  and  approaches  a  limiting  value  as  t  -0.  Whalen^^ 

8 

reported  a  similar  result  for  the  displacement  thickness  for  Pr  =  uu  =  1. 

It  is  very  difficult  to  obtain  exact  numerical  solutions  of  the 
laminar -boundary- layer  equations  for  values  of  3  grea.er  than  2.  The 
reason  is  simply  that  the  singular  behavior  of  f  /(T|)  near  the  wall  as 
P  -•  ®  becomes  dominant  even  for  moderate  values  of  p.  Conversely,  the 
results  obtained  for  P  -  ®  should  be  good  representations  of  the  be¬ 
havior  of  the  boundary  layer  for  large  but  finite  values  of  P.  One 
of  the  imp or tant  results  we  wish  to  demonstrate  is  chat,  by  combining 
the  exact  numerical  results  for  P  5  5  given  in  Tables  1  to  5  and  the 
limiting  solutions  for  P  -*  ®  given  in  Tables  6  and  7,  it  is  possible 
to  estimate  accurately  the  skin-Irictio.i  and  heat-transfer  derivatives 

t" ,  g',  and  0;  for  all  positive  values  of  p. 
w  w  w 

The  skin-friction  results  for  P  *  ®  are  displayed  in  Fig.  10. 

The  influences  of  the  two  parameters  w  and  a  on  the  inner  solutions 
for  8^(0)  are  seen  to  be  relatively  small.  For  0.5  s  id  £  L0  and 
0  5  a  5  1.0,  the  value  of  aj(0)  may  be  found  by  interpolation  co 
better  than  0.25  percent. 

The  skin- friction  parameter  £_"()  *^(tw/t  )  approaches  the 

limit  of  a,(0)  as  P  -*  ®.  The  difficulties  that  were  observed  previ- 
o 

ously  in  calculating  exact  solutions  for  p  >  2  imply  that  this  limit 
is  approached  very  rapidly  with  increasing  P.  This  supposition  is 
borne  out  by  Figs.  11  to  13,  wheie  the  skin  friction  parameter  is 
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shown  as  a  function  of  0  .  The  limit  parameter  0  Is  suggested 

by  the  ordering  procedure  used  to  obtain  the  Inner  and  outer  equations. 
Solid  lines  Indicate  exsct  numerical  solutions  and  dashed  lines  are 
extrapolations. 

Figure  11  Illustrates  the  approach  of  the  skln-frlctlon  parameter 
to  Its  asymptotic  limit  s^(0)  for  different  wall  temperatures.  The 
limiting  value  is  approached  most  rapidly  for  high  well  temperatures. 
This  result  is  to  be  expected  from  the  behavior  of  the  outer  equations. 
Large  values  of  t  Increase  the  magnitude  of  the  velocity  difference 
across  the  inner  layer,  wherea«  for  t^  -*  0,  the  distinction  between 
the  inner  and  outer  layers  breaks  down  and  no  proper  limit  la  obtained. 
It  has  also  been  found  empirically  that  exact  solutions  are  more  diffi¬ 
cult  to  obtain  for  large  t^. 

The  approach  of  the  skin-friction  parameter  to  its  limiting  value 
with  1  '.creasing  0  is  illustrated  in  Fig.  12  for  several  values  of  the 

sweep  parameter  t  .  From  a  numerical  point  of  view,  the  accuracy  of 

0 

the  present  extrapolation  procedure  increases  with  increasing  sweep. 

It  is  of  interest  to  examine  the  behavior  of  the  inner  and  outer 
equations  with  mass  inject' an  at  the  wall.  Applications  of  these  re¬ 
sults  include  mass-transfer  cooling  of  rocket  nozzles  (Back  and 
Witte)^)  and  blunt  hypervelocity  vehicles.  The  outer  equations  de¬ 
termine  the  heat- transfer  derivative  9j,  and  they  contain  the  boundary 
condition  f  (T1  -*  0)  ■  f  .  Thus,  surface  mass  transfer  (f  <  0)  acts 
to  reduce  and,  consequently,  surface-heat  transfer  even  in  the  limit 
0  -*  *.  On  the  contrary,  the  skin  friction  derivative  f^  is  found  from 
the  inner  equation  for  ?J(0),  Eq.  (147),  a  the  solution  of  this  ••  ua- 
tion  is  independent  of  the  value  of  f^.  In  highly  accelerated  flows, 
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therefore,  the  effects  of  blov  Lng  on  skin  friction  become  negligible 

in  the  limit  3  -•  <*>  with  f  fixed. 

w 

Figure  13  illustrates  the  behavior  of  the  skin-friction  parameter 

as  a  function  of  the  injection  parameter  f^.  The  limit  is  approached 

-1/2 

smoothly  with  decreasing  values  of  (3  for  all  f  ,  and  accurate  esti¬ 
mates  of  f*  may  be  obtained  for  all  {3  by  comparing  the  exact  solutions 
for  j3  5  5  and  the  inner  limit  solutions  of  Fig.  10. 

The  wall  heat  transfer  is  related  by  the  modified  Stewartson  and 
riowarth-Dorodnitsyn  transformations  (Eqs.  (10)  a-  1  (11))  to  0j.  Inas¬ 
much  as  Sq  anf  g  represent  the  complete  solutions  tor  9  and  g  as 
9  the  values  of  9^(0)  and  g^(0)  represent  the  asymptotic  limits 

of  0'  and  r  7  as  (3  —  ®.  .  t  should  be  emphasized  that  9  7  and  g7  be- 

w  °w  w  °w 

come  independent  of  (3  as  3  -*  °°,  demonstrating  that  the  heat  transfer 
predicted  by  a  local  similarity  analysis  ip  highly  accelerated  flows 
approached  a  limiting  value. 

Figure  14  shows  the  typical  behavior  of  the  heat-transfer  deriv- 

tive  97  with  increasing  values  of  9*  The  approach  of  9 7  to  its  limit- 
w  w 

ing  value  for  9  -*  ®,  is  seen  to  be  smooth  (at  least  for  the  case 

(4) 

of  ct,  =  0)  monotonic.  In  an  earlier  paper,  7  we  demonstrated  that 

J. 

the  proper  parameter  to  use  in  comparing  different  heat-transfer  cal¬ 
culations  is  97[(1  -  t  )/(t  -  t  )].  For  c.  =  0  and  t  =  1,  the 

wUN  w  aw  w  J  1  s 

adiabatic  wall  temperature  t  is  unity  for  all  Pr  so  that  the  heat 
transfer  parameter  reduces  to  9j.  Although  we  have  nc  been  able  to 
prove  it  analytically,  it  appears  that  t  =  1.0  for  all  values  of 
Pr,  ts>  id,  and  a  in  the  limit  9  =  <*>•  This  is  a  very  s'-^prising  result 
and  should  be  examined  further. 
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The  boundary- layer  displacement  thickness  6  is  defined  by  the 
relation  [Eq.  (79)1 


r”  /  P“  V  /T  \ 

/Tv 

J  (*  -  )iy  -  j  L ) 

‘l  -If  1 2 

o'  p  u  /  p  u  r,  ' T  ' 

e  e  e  e  k  e 

'T  f 

L  0 

where  the  integrals  1^  and  are  given  in  the  list  of  symbols.  For 
P  »  ®,  the  integral  1^  is  identically  zero  and  the  integral  [Eq. 
(75)1 


I 


2 


«® 

f'(l  -  f')  d"H 

O 


is  recorded  in  Table  6.  In  the  absence  of  sweep,  the  velocity  profile 

is  monotonic  and  0  <  f '  <  1,  so  that  the  integral  is  positive  and 

the  displacement  thickness  is  negative.  With  sweep  (t  <  1) ,  there 

s 

is  often  an  overshoot  in  the  velocity  profile,  so  that  f '  >  1  for  some 
range  of  T),  and  I  becomes  negative.  With  sweep,  therefore,  the  dis¬ 
placement  thickness  may  be  either  positive  or  negative,  depending  upon 
the  particular  parameters  being  considered.  Numerical  comparison  be¬ 
tween  the  calculated  results  for  moderate  P  and  the  present  results 

* 

for  p  =  <»  suggests  that  6  raonotonically  decreases  with  increasing  p. 

In  the  special  case  when  P  =  ®  and  t  ■-  t  =  1,  f '  =  1  and  h  is  of 

w  s 

-1/2 

the  order  of  P  times  the  scale  of  the  boundary  displacement  layer 
thickness  for  P  =  0. 

The  numerical  results  giv-m  .n  Table  5  display  the  variation  of 
boundary- layer  properties  for  various  Prandtl  numbers  near  unity. 

These  results  may  be  used  to  estimate  recovery  factors  for  Pr  >  1  by 
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combining  Che  present  results  with  the  asymptotic  solution  of  Narasimha 
and  Vasantha^^  for  Pr  »  1.  Although  these  authors  explicitly  solved 
the  problem  of  a  flat  plate  (P  *  0)  with  to  ■  1  and  low  Mach  numbers, 
there  is  no  reason  why  a  similar  analysis  could  not  b<i  conducted  for 
P  >  0  and  general  compressible  flow.  As  Narasimha  and  Vasantha  demon¬ 
strate,  interpolations  accurate  to  better  than  3  percent  in  the  re¬ 
covery  factor  can  be  made  between  first-order  asymptotic  solutions  for 
Pr  »  1  and  exact  numerical  results  for  Pr  of  order  unity. 

In  concluding  this  discussion,  it  is  advisable  to  point  out  that 
the  numerical  values  listed  for  the  wall  derivatives  f",  9/,  and  g,' 

W  W 

allow  reconstruction  of  the  complete  velocity  and  total  enthalpy  pro¬ 
files  by  standard  numerical  integration  techniques.  Whereas  the  bound¬ 
ary-layer  equations  themselves,  including  the  boundary  conditions  at 
T|  ■  0  and  T|  -*  <*>,  represent  a  two-point  boundary-value  problem,  one  may 
solve  the  equations  as  an  ititial-value  problem  if  the  wall  deriva¬ 
tives  £*,  0  and  g/  are  known.  Furthermore,  by  using  the  wall  derl- 

W  W  vf 

vatives  given  here,  multi-term  expansions  of  the  boundary- layer  equa¬ 
tions  may  be  constructed  for  small  11.  This  property  is  useful  in 
cases  where  additional  "nonclassical"  behavior  occurs  near  the  surface, 
for  example,  in  MKD  boundary  layers  where  an  electrostatic  potential 
sheath  exists  for  small  T). 


Fig. 6— Variation  of  the  Reynolds  analogy  function  with  pressure  gradient 


Fig. 7  Variation  of  the  Reynolds  analogy  function  with  pressure  gradient 


{l  +[g  " (0)  sin  A/ f  "  (0)  cos  A]2} 


6 


in  friction  with  injection  for  large  0 
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Fig.  14— Variation  of  heat  trcnsfer  with  wall  temperature  for  large  0 
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